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Abstract 

We demonstrate that Aharonov- Albert- Vaidman (AAV) weak values have a direct relationship with the 
response function of a system, and have a much wider range of applicability in both the classical and 
quantum domains than previously thought. Using this idea, we have built an optical system, based on a 
birefringent photonic crystal, with an infinite number of weak values. In this system, the propagation speed 
of a polarized light pulse displays both superluminal and slow light behavior with a sharp transition between 
the two regimes. We show that this system's response possesses two-dimensional, vortex-antivortex phase 
singularities. Important consequences for optical signal processing are discussed. 
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The role of measurement in quantum mechanics has posed challenges for physicists for nearly 
a century. Standard quantum mechanics represents ideal measurements as the action of a Hilbert- 
space operator on a state. In an ideal or "strong" measurement, the operator completely collapses 
the state into an eigenstate, with an expectation value within the range of the operator's eigenval- 
ues. Strong measurements yield maximal information about the state of the system and correspond 
to most laboratory measurements. 

Contrary to traditional quantum measurement theory, certain experiments can return results that 
are far beyond the range of an observable's eigenvalues |l|]. These so-called "weak" measurements 
barely perturb the state, and have been less studied in quantum mechanics because they yield very 
little information. In their seminal work, Aharonov, Albert and Vaidman (AAV) showed that para- 
doxical results can be obtained if one preselects an initial state, performs a weak measurement, 
and subsequently a standard strong measurement which postselects a particular final state. Al- 
though the action of the weak measurement has little effect on the global state of the system, it can 
have a significant impact on one or more of its components. Weak values have been theoretically 
discussed for spin- 1/2 particles and experimentally demonstrated in optical systems involving the 
angular deflections of polarized light beams passing through birefringent prisms III 1^1]. 

Recently, using the approach of quantum-trajectory theory, it has been demonstrated that post- 
selected weak values can be thought of as quantum correlation functions, which can be used in 
quantum optics and related areas [5]. In this paper, we instead start from the original derivation of 
the AAV effect and extend it further. 

We show that weak values appear in a natural way in a much broader range of situations. 
We demonstrate that a system's response function (or scattering matrix) has a direct relationship 
with AAV weak values. As an illustration, we present an experimental study of a frequency- 
dependent, polarization-sensitive, optical system based on a two-dimensional (2D), birefringent 
photonic crystal. The response of this system can be described by an infinite number of unbounded 
weak values which physically correspond to measurements of the polarization value or spin of the 
photon along different axes. Unlike previously studied weak measurements, this system maps 
polarization states onto collinear spatial position along the direction of propagation, resulting in 
unbounded values for the time-of-flight of a polarized optical pulse. The system displays both 
superluminal and slow light [8] behavior with a sharp transition between the two regimes 
and gives us total control over the speed of light. It may lead to several important applications 
in optical signal processing J^], such as optical delay lines with polarization-tuned characteristics 
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and dynamical optical clock multipliers. This new approach to quantum weak values should also 
open the door for many new precision measurements of very tiny quantities JilQ]. 

The complex response of a linear system parametrized by two general variables (p, 77) is given 

by 

f(p,ri) = (xl> f \U(p,ri)\4> in ) (1) 

where U (p, 77) is a unitary operator describing the evolution of the system, and \ipi n ) and \ipf) 
are normalized input and output states. To apply the AAV formalism, we must express Eq. [l] 
in terms of the preselection of an initial state, a weak measurement, and a strong measurement 
or projection onto a particular final state. The difficulty lies in the fact that U (p, rj) cannot be 
generally classified as a weak operator because it may have a strong impact on \tpi n ). However, 
we are free to decompose U (p, rj) into a large or "strong" part and a small or "weak" part. To 
do so, we Taylor expand the general operator U (p, 77) around an arbitrary point (p , r] ) for small 
deviations Sp and Sr], and separate U (p, 77) into a weak operator C/ w (Sp, 5rj; p , 770) and a strong 
operator Us = U (p ,rj ) obtaining 

f (p + Sp, r/ + 8ri) « 

(^/ \4 ) in(po,Vo)){l+'i[A p } w 5p + i[A v } w 5r]} (2) 
where ijj in (p , 770)) = U s \ip in ), and [A p ] w and [A v ] w are given by: 



(i>f\Aj(p 




4>in (f>0,Vo)) 




4>in (Po,Vo)) 



[As (Po, Vo)] w = — ^ (3) 

with Aj (p , 770) = -i [djUs] U^ for j -> p, 77. 

The strong operator Us can be viewed as a preselector which generates the state ip in (p , 770) 
from \ipin), and U-w (Sp, Si]; p , 770) is a weak measurement operator which acts on this preselected 
state. Eq. [3] then implies that [Ap]^ and [^4r?] w are AAV weak values associated with the oper- 
ators A p (p , 770) and A v (p , 770). These weak values [A p ] w and [A, n ] w can be unbounded if the 
corresponding transfer function has a singularity, i.e., a point for which [ipf ^) in (p , r] )^ = 0. 

This analysis is valid for classical and quantum-mechanical systems. In the quantum mechan- 
ical case, we must reinterpret 5p and 5i] as operators rather than deviation parameters, while, the 
preselected state is a superposition of eigenstates, sharply peaked about p an d 770 such that the 
expectation values of Sp and Si], appearing in Eq. El are small relative to po and 7/0, respectively. 
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FIG. 2: Experimentally measured phase delay vs. frequency for microwaves transmitted through an 18-layer 
photonic crystal positioned at different angles. 

Experimentally observable quantities are inferred from an indicator or pointer on a measuring 
device. In the general case, 5p and 5r) are canonical variables describing two independent measur- 
ing devices. The pointers of these measuring devices are the conjugate momenta to these canonical 
variables, and the interaction weak operator XJ W couples them to the physical observables A p and 
A v . As a result, the measurement process shifts the pointer by an amount equal to the weak value. 
Our analysis shows that there is always a well-defined connection between measurement pointers 
for weak values and the response function of a system. To make this explicit, we invert Eq. El 
obtaining 

[Aj {po,m)) w = 

dj [argf (p , Vo)\ ~ idj [in f (p , Vo)\] ■ (4) 

In words, the local gradients of the phase and the logarithm of the magnitude of a response function 
are measurement pointers for the real and imaginary parts of weak values. 

To illustrate this formalism we apply it to experimental measurements of the transmission of 
normally incident, linearly polarized microwaves through a photonic crystal. These crystals, con- 
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FIG. 3: Experimentally measured group delay vs. crystal rotation at uj = 16.7 GHz (a) and pointer for 
vs. frequency at (3 = 7r/4 (b). 

structed using a previously published method lioll . are highly birefringent Jill [l^i . Our crystal 
can be rotated in the a;2;-plane by an angle (3 with respect to the ^-axis and its fast axis is parallel to 
x at (3 — 0. Microwaves polarized along z are generated by a HP 8720A vector network analyzer 
(VNA) and detected by a polarization-sensitive horn (see Fig. 1). The two degrees of freedom oj 
and [3 allow the transmitted polarization to lie anywhere on the 2D Hilbert space of polarizations 
(Poincare sphere). If the crystal is positioned at (3 = it/ A, there is no detected signal at uj s , the 
frequency for which the crystal behaves as a half waveplate, because the transmitted EM wave is 
entirely x-polarized. 

In Fig. 2, we display experimental results for the phase delay as a function of frequency for 
microwaves transmitted through an 18-layer photonic crystal (u s = 16.7 GHz) positioned at vari- 
ous angles. The phase delay exhibits strong dispersion near (3 = 7r/4. Surprisingly, for P < it/A 
the dispersion is anomalous, whereas it is normal for (3 > 7r/4. Some insight into this behavior 
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FIG. 4: Theoretical contour plot of arg T(u>, (3) showing the singularities and the surrounding phase. The 
arcs connecting the singularities correspond to 2ir phase slips. 



can be found from a general topological argument. If the parameters (u, (5) are swept to trace out 
a closed loop around the singularity (zero-transmission point), the phase can change by an integer 
multiple of 2n over a complete cycle. The phase may therefore show extremely rapid and even 
discontinuous changes as these parameters are adjusted through paths that pass close to or over the 
singularity. 

To apply our framework to the experiment, we construct the unitary operator U (u, (3) in terms 
of the usual Pauli matrices a k , finding 

U(u, 0) = e xp[i0 TE (i4p) + z0 TM (l^)] (5) 

where 0te, 0tm are the birefringent phases for TE and TM polarizations and k — > 1,2,3. The in- 
dependent vertical and horizontal polarization states 1 1 ) and 1 2) are parallel to z and x, respectively, 
and cr 3 |1) = — |1) and cr 3 |2) = |2). U(u, (3) is found by a simple rotation, and the corresponding 
weak operators follow naturally from the Taylor expansion. 

Our weak values are associated with polarization measurements of the incident radiation. For 
our experimental system, the imaginary part of Eq. — > oS) vanishes for points along the line 
(u> s ,Po). Thus along the line (a> s , (3 ) the pointer for [A^]^ is simply the frequency derivative 
of the transmitted phase, or the group delay of a pulse with narrow spectral width Su and center 
frequency lu s . The weak value [A u ] w is coupled to a spatial position pointer and manifested as a 
time- of- flight. A calculation of the weak value yields 

[A^Us, f3 )] w = sec(2/3 o )[cos 2 (/3 o )<9 a ,0TM(w s ) 
6 



-sin 2 (A))<9u0TE(u; s )]- (6) 

By taking numerical derivatives of the measured phase delay with respect to the frequency, we 
generate experimental data for the pointer for [A^ (u a , Po))w m ^S- ^(a) we show the results at uj s 
as a function of the waveplate angle. For comparison, we also show theoretical curves generated 
using Eq.|6]and independent measurements of 4>tm (oj) and (f) TE (co). 

The group delay measurement pointer shows both "fast" and "slow" light behavior; on one side 
of the singularity at /3 = 7r/4, arbitrarily large, negative group velocities are possible {(3 < vr/4), 
while on the other side {(3 > 7r/4), pulses can propagate with arbitrarily small group velocities. 
The pointer takes different positive, non-degenerate values at ft — and (3 = tt/2. As a result, this 
passive system can be used as a compact optical delay line whose properties can be tuned by small 
adjustments of the waveplate angle and the birefringence. Since the birefringent characteristics 
of a photonic crystal can be tailored by changing its geometry, air-filling fraction, and material 
composition [12], this type of delay line can be adapted to suit a wide variety of applications. 

The weak value [Ap (uq,(3 s )] w is also a polarization value or spin of the photon. Near the 
singularity, the complicated form of the operator Ap (u; ,/3 s ) reduces to 2<r 2 = L y /h to leading 
order in cu . Thus, in this vicinity, the weak value is simply the total helicity of the light about the 
y-axis. Along the line (u> , (3 a ) where (3 S = tt/4, the pointer for the weak value reduces to the real 
quantity dp arg T (uq, (3 s ) . An explicit calculation of the pointer on this line produces the simple 
form 

[Ap (w , p s )] w = 2 tan [</>_ (cu )} . (7) 

If 0_ (o) ) — \ [4>te{u) — 4>tm(u)] = or 7r, Eq. |7] yields the degenerate eigenvalues 
[Ap (ujQ,P s )) w = 0; however, the pointer assumes arbitrarily large positive and negative values 
near the tangent singularity, where it corresponds to a measurement of the total helicity of the 
light about the y-axis. Since the number of photons in light is of course finite, this unbounded 
weak value would imply that the basic unit of angular momentum carried by each photon can be 
anything. In Fig. 3(b), we display experimental data for this pointer as a function of frequency 
with (3 = 7r/4. 

One may wonder why the crystal itself cannot be used as a measurement pointer for this weak 
value since the angular momentum of the crystal about the y-axis is conjugate to (3. As a classical 
object, it cannot exist in a superposition of angular states. Therefore, its angular momentum is 
uncertain on the quantum scale and will not reflect the weak value. This point emphasizes the 
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value of Eq. local derivatives of the response function can always be used as measurement 
pointers for weak values, regardless of whether the experiment is classical or quantum-mechanical. 
Similar reasoning applies to the group delay pointer for [i^J^,, since a pulse is simply a classical 
superposition of different frequencies: if the weak measurement is performed with continuous- 
wave light, the time-of-flight is not a well-defined quantity, but can be directly inferred from 
W argT (u ,Po) . 

In Fig. 4, we display a theoretical contour plot of argT (to, (3) generated using simple lin- 
ear models for cf) TM (u) and 4>te (^)- Singularities are expected at all points (uj n ,f3 m ) where 
0_ (uj n ) = (2n + 1) | and (3 m = (2m + 1) | for integers n and m. Using these models, we dis- 
play eight singularities, and the topology of the surrounding phase. The real parts of the weak 
values [A w (uo 0) (3o)] w and [Ap (uq, (3o)] w at any point are given by the local gradient of the phase. 
Between the singularities, there are saddle points, i.e., points where the weak values are least 
sensitive to deviations in the parameters. The global structure is that of a vortex- anti vortex lat- 
tice, with adjacent singularities possessing opposite topological charges. This structure implies 
that the system is relatively insensitive to perturbations in the operator U (to, (3) because the spon- 
taneous appearance or disappearance of a singularity would have a discontinuous effect on the 
global phase. Thus, the topology must possess an equal number of vortices and antivortices; as 
the unitary operator is continuously perturbed, vortices and antivortices can either approach each 
other and annihilate or separate to infinity 11311 . Finally we note that while [AJw an d [A?W WQ 
the natural weak values of this system, arbitrary superpositions of the two weak operators lead to 
infinite number of other weak values whose physical interpretation is less clear. 

In conclusion, we have shown that response functions can be treated with the AAV formalism, 
although the physical interpretation is specific to the particular physical system under study. If 
the response possesses singularities, one obtains weak values which can be far outside the ranges 
spanned by the eigenvalues of the corresponding operators. We have demonstrated this using 
an optical system exhibiting singular points in a two-dimensional state space characterized by 
two independent weak values and an infinite number of superpositions thereof. This system has 
significant implications for optical communication applications requiring compact, tunable delay 
lines and clock multipliers. Although our experiments were performed using classical signals, the 
results apply to quantum-mechanical experiments as well. 
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